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ABSTRACT 



Geostron’iic relative vorticit;/ is confuted along three latitude 
circles generally encompassing the United States, using a Fourier analy- 
sis of the 500-rab height field for the period 3-6 November 19^1 . 

A statistical comparison of vorticity obtained in this manner with 
it*s square grid computation is made. 

The writers are deeply indebted to Professor F. L, Martin of the 
United States Naval Postgraduate School for his suggestion of the topic 
and his continued help throughout the investigation and during the 
preparation of this paper* 
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Symbol Definition 



% 


vertical component of geostrophic relative vorticity 


g 


acceleration due to gravity 


f 


coriolis parameter 


z 


contou.r height 


a 


radius of earth 


0 


angular velocity of earth 
latitude 


a r 


difference in latitude 


* 


longitude 




4 JGL - 1; f/ i 


'IT 


zonallv averaged contour height 


Z " 


perturbation of zonal contour height 


z 


space -mean contour height 


d 

H 


grid mesh = 381 km 
wave number 

latitude circles numbered 1 through 27 grid lengths from the 
north pole 


C n 


amplitude of nth wave 


©n 


phase angle of nth wave 


m 


map factor 


r 


radius of earth at latitude <b , ^ - d Cj; ‘ 


J 

F/V//F 


vertical component of geostrophic relative vorticity as com- 
puted by Fleet Numerical Weather Facility 


F 


F test for equality of variances 


p 


number of predictors (independent variables ) 


u g 


tro r nic;-liy -computed zonal component cf the wind 


r 


i— o.l wl.. 1 v r. . t .ed geostrophically 
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Symbol 



Definition 



RKS 

R 

'n 



root mean square 

multiple correlation coefficient 
number of cases in sample 
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1 . Introduction. 



The advent of the electronic computer coupled with increased data 
coverage has removed the primary obstacles to the use of harmonic analy- 
sis as a useful instrument in atmospheric research . 

Use is made of a Fourier analysis of the 500-mb height field. A 
computer program devised by G. Arnason and M. Reese, and made available 
by Fleet Numerical Weather Facility"' was used to obtain the harmonic 
analysis around latitude circles at distances of 1 through 27 grid lengths 
d = 381 km from the pole. An equation for relative vorticity using the 
results of the spherical harmonic analysis is derived and tested statisti- 
cally against the more usual numerical computation of geoctropbic relative 
vorticity made using a square grid. 

The problem undertaken was to compute relative vorticity at 500 mb 
using the peostrophic vorticity equation in its spherical coordinate form 
[1, P- 358], 
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In equation (l ), z was replaced by its expression in terms of a Fourier 
analysis 

7 ~ i t £[_C,, -e\)'J - f A- /* (2) 

n ; ( 

Here z is the zonally averaged contour height, whereas z*‘ is the perturba- 
tion quantity which remains when the z field is subtracted from the actual 
height field. The computer printout included maps of each of these three 
fields at 500 mb for the period under investigation, 3 through 6 November 
1961. Geostrophic relative vorticity was computed by the above mentioned 
technique for 17 longitudes along each of three latitude circles, generally 
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Henceforth referred to as FNVJF 
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encompassing the north-south extent of the United States. This gave a 
total snmnle size of 204 computed relative vorticities which were then 
correlated with square-grid computations of relative vorticities obtained 
from FNUF charts of absolute vorticity for exactly the same set cf points 
and times. 

2. The geostrophic relative vorticity in terms of a harmonically analyzed 
contour field. 

The adequacy of the Fourier representation of meteorological para- 
meters will not be discussed here; the interested reader is referred to 

L2J. 



It was previously mentioned that use was made of the FNWF computer 
program for a Fourier analysis of the 500 mb contour fiold. A brief de- 
scription of the theory and numerical evaluation of the Fourier coeffi- 
cients is given in [3 j and only the applicable portions of the program 
will be mentioned here. 

The zonal mean height z is defined as 
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(3) 



and its zonal perturbation, in accordance with equation (2), is 

2 -I (4) 

Expanding z in a Fourier series one obtains at a given latitude 

•i 

2 ~ 2- ^ 11 [A. 1 l3 n C ' c 5 0 (5) 

n - 1 

where A and B are the coefficients of the nth harmonic functions in the • 

Yi r\ — 

Fourier expansion of z, and M is the total number of harmonic waves neces- 
sary to resolve the contour field. Replacing and in (5) by an am- 
plitude factor C ^ and a phase angle 0 , (5) may be revrritten as 



2 



'? 



( 6 ) 



i 



l 



c 

n 






/// t t’ A — t: M ) 



w*>ere 



c„©a>v 



ana fan o„ -B„/ a„ 



Equation (6) has been programmed for the CDC-1604 digital computer by 
FNUF, and the program includes: a contour grid map of the zonal pertur- 

bation field z*, and a tabular list of the Fourier parameters C n , © n and 
z for each of 27 latitude circles. The last three parameters are actually 
computed in relation to a zonal band of width equal to the grid mesh d of 
the FNUF map grid, but may be considered applicable at the center of the 
band . 

Using spherical coordinates h and (b (longitude and latitude, respec- 
tively) the derivative operators 
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apply, and we have 
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4/r : ~ Cnt <P 



(7) 



If nov/ z is replaced in (7) by its expansion in a Fourier series as given 
in equation (2) we obtain 

h ~ of' l c \_L. 1 (1 111 ( ' " 0, =Y 
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+~ [_L jlh k i,A ~ ■*" u i < ~* c ; f (s) 

Remembering that C ^ and p are independent of longitude 4 but not of 
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latitude .f , we may perform the indicated differentiations and combine 
terms, with the result that j f in spherical coordinates becomes 

/, - £, [-.-.j tQsmlr.A *«.)( y" - -,(]'■■■) )j 
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Equation (9) may be shown to have the form 
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equation (l( ) may be approximated by 
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Here the superscript bar denotes an average with respect to n over all 
significant waves n = 1 , , M. 

The FNWF maps are stereographic naps true at latitude 60N, and else- 
where the grid mesh, d = 3&1 km, is related to the true distance d T by the 
relationship d = md_, where m - (l + sin 60)/(l + sin,|>) and is the so- 
called map factor. 

The nartial derivatives in (ll ) nay be evaluated by centered finite 
differences. The numerical values of C n and v;ere simply token from 
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the tabular list of computer output of the harmonic analysis at the vari- 
ous latitude circles. For example, d / i ' may be written using 
finite differences between adjacent latitude circles (for which a a cj-, 1 An ) 
in the form 
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The latitude circles of the Fourier analysis employed in the compita- 
tion of J<j corresponded to N = 12, 15, and 18 grid mesh distances from the 
pole, which in turn corresponded to the latitude circles: 48. 9H, 38. 6N, 

and 28, 3N (see fig, l). However, in the computation of first and second 
derivatives of the nth Fourier wave parameters, these parameters at lati- 
tude circles I'M and N+1 are also needed. 

The term vL/ f\ is more tractible in the present case if the finite- 
difference substitution ri A - 2 . d/m is made; thus, the centered deriva- 
tive A in finite differences, assumes the form 




where T - c/ c--: j is the radius of the latitude circle, and the subscript 
i+1 is one grid mesh eastward of the ith point alopg the latitude circle. 

Including the form of mV.;A. given in (12), equation (ll ) may now 
be written in the more compact form 
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Fig. 1. Grid overlay depicting latitude circles N=12, 15, IS and the grid mesh length, d. Geostrophic relative 
vorticity is con put ed at the intersections of the latitude circles with all 5 'longitude lines lying 
within the area bounded by the double "line meridians. 



The averaging with respect to n in the expressions for W and is 
assumed to be valid, according to the mean value theorem. Actually how- 
ever the average was approximated by a simple arithmetic mean of the con- 
tributions for n = 1 , . M. Thus it was necessary to compute contri- 
butions for each n and divide the result by the number of waves. 

Van Mieghem |^2 j points out that waves of space number higher than 12 
to 17 contribute little to the variance of the contour-height field. Con- 
sequently, the Fourier analysis was truncated at M = 12 at the two higher 
latitude circles, 1 although wave numbers 1 through 15 were employed at the 
lowest latitude. 

The coefficients and W 2 are constants for a particular time and 
latitude circle, and can be computed directly from the tabulated harmonic 
data. These values are listed in table 1 . However, in applying equation 
(13) it was necessary to interpolate values of z"" from the pert”.rbation 
map at the points indicated in fig, 1 along the specified latitude circles 
These points were arbitrarily chosen at integral multiples of five degrees 
of longitude, encompassing the west-east extent of the United States, 
Therefore, for ease in selecting data points, a grid overlay was con- 
structed presenting all latitude circles N?=11 , » . , ,,19, drawn so as to in- 
corporate the map factor. 



At the two higher latitudes, H, was negative ensuring that lows 
corresponded to positive relative vorticity. At the lowest latitude 
(N= 18 ), it was necessary to use M=1 5 in order to ensure that this condi- 
tion was met. 
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AID values 



TAB IE 1 : Values of U, and W x to be applied in equation ( , . 

. 1 _ i 

are in units of 10 sec ft . 



1“ "" LATITUDE CIRCIE 

I 1 



DATE 


1 n 

N = 12 


1 


3 

II 

-1 1 

V/l j 

1 

1 


r 


N = 18 




W = -.1467 


W 


= -.1635 


W 


= -*2696 


3 Nov. 1961 














w = -,0079 

1 


• u 

1 


= -.0370 


w 


= -.0303 




W = -*1490 


w 


= -.0617 


w" 


= -,0420 


4 Nov. 1961 


i 












w = -.0311 


w 


= .0316 


w 


= .0138 




w = -,1355 


w 


= -*1013 


w 


= -.1819 


5 Nov. 1961 




11 = -.0398 


w 


= -.0230 


• W 


= .1682 




w = -,6551 


w 


1 

"i 

1 ! 

* 

£ 

O! 

°! 

1 

I 


u 


= -.3510 


6 Nov % 1961 


? 


w = .0757 

L___ _ 


w 


= .0349 


w 


= -.7610 



Mote that the term ( 2 u.., - t V-- in both (ll ) and (l 3 ) also varies 

from point to point around any latitude circle* However, the variation of 
this term was so small that the mean zonal wind for the latitude was used . 
The mean zonal ’wind was computed using the latitudinal gradient of 
zonal mean height z, that is from 

~ ^ ^ ^ ' * / t * _ 1 rfjJ. -j 

3 r zd ’ ; 

again using finite differences* 

Along each of the three latitude circles, 17 values of W, z* and 
W (z * ( ) were obtained for each of the four days, giving a total of 
204 separate values of these quantities, 

3. Statistical results. 

For precisely the same set of ooints and days as described above, 
corresponding values of */ere obtained from grid maps of geostrophic 
relative vorticity made available by FNTJF . Their maps essentially pre- 
sented the field of J j computed using their square-grid mesh, according 
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to the usual formula 
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(14) 



where z is the four-point s pace - mean contour height and d = 381 km at 
latitude 6 ON* 

Since the averaging process in equation (ll ) is somewhat approximate, 
it was felt that better verification would result if the relationship be- 
tween J r and U z* and W„ ( z ?' -z f , ) was written 

f - rv *y r t l a, a / L i 



J, 
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A. t- 



)1 ( 15 ) 



Here A^ , A. , and A g are regression coefficients^ to be determined by a 
least squares technique applied to the 204 simultaneous values of W, zf 
and W- (z* ( -z5' ( ). Because of occasional large gradients of '« rNVVf over 5° 
longitude intervals, the values of / were smoothed by the three-point 
binomial running mean formula 




* 'h + *<>,) /"/ 



where is any value of j ( interpolated from the FNUF vorticity-grid 
map at each point. The regression equation actually involves the correla- 
tion of the smoothed J ^ .. values with the independent variables already 
noted . 

Using the BIMD 06 statistical program as adapted to the CDC-1604 
computer, it was necessary to type 204 IBM cards containing the smoothed 
J ^ values and the simultaneous values of the independent variables of 
the right side of equation ( 15 )- Best-fitting values of the regression 
coefficients were 



These regression coefficients are non-dimensional since the units 
of U. zf and VI (z* -zq. ) are in sec-” 1 . 
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A c = -8.1719, 



A 1 = 0.4585, 



and A ? = 0.0624 



The standard error of esti.mate"' and the multiple linear correlation co- 
efficient for the regression equation were 20.3523 and 0.7632 respectively. 

Panofsky and Brier [4, p. 113] present an F test for the analysis of 
the variance explained by the regression equation. The F value may be de- 
fined by the ratio of the explained variance to the residual variance 
relative to the regression equation. This formula is 



ft* 

a rv*) p 



(16) 



With p = 2, the number of independent variables, and 6} = 204, the number 
of random samples, the value o.f F from equation (l6) turns out to be 
139.33. However , assuming normally distributed independent and dependent 
variables, the critical F at the 0.99 significance level with 2 and 201 
degrees of freedom is 4.71. Thus the regression equation (15) is signi- 
ficant at a very high confidence level. 

The partial correlation coefficients of U zr and W„ (z." - 2 * ) were 
0.7620 and 0.0382, respectively, indicating very little contribution to 
the regression equation for the latter variable. An intuitive explana- 
tion for the insignificant contribution of W^ (z* -z^ ) may be obtained 
from fig. 2, The differences, z?' : '. -z?f , taken over the pressure centers 
would be zero, whereas in actuality the largest absolute values of J ^ 
would be found at these points. likewise for a sinusoidal contour wave 
pattern maximum magnitudes of 2 ‘ :: ' - 2 * would exist near contour inflection 

4*i A. - * * 

points, where i } - . 

The possibility of "automatic" correlation, viz, correlation between 
the two independent variables, was examined and none was found to exist. 



All values of are in units x 10“^ sec -- ' 
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Fig. 2. Distribution of along a latitude circle between contour maxi- 
mum and minimum, in the case of a sinusoidal contour wave. 

The individual latitude circles were statistically examined in a like 
manner. The two northernmost latitude circles gave a multiple linear cor- 
relation coefficient of 0.84, whereas there is no significant coi'relation 
between J and the two independent variables W.zf and W, (z* -z* ) at 
latitude 28 ,3N. This inconsistency at low latitudes is attributable to 
"noise" in the Fourier representation of the contour field. The "noise" 
is due to the small variance of the height field at this latitude, coupled 
with the relatively large RMS error possible in the 500 mb analysis in the 
subtropics during the period of study. The use of second derivatives of 
the Fourier parameters also increases the amount of "noise" generated in 
the determination of 1 , . 

4. Conclusions. 

Theoretically there should have been a large positive multiple linear 
correlation between j and W z V and W. (z;* -z. ,), This proved to be the 

r i a. i »-/ x 

case at latitudes 48. 9N and 3S.6N where the contour field was well resolved 
by twelve harmonic waves. Fowever, at the lowest latitude circle, 28. 3F, 
the several sources of errors listed at the end of section 3 could, and 
probably did, influence the low correlations that were obtained. These 
were basically due to the lack of resolution of the contour field by the 
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Fourier analysis. However, the cane reasons would cause the square-grid 
commutation of.) > to be subject to large errors. 

Based upon the results at the two higher latitudes, it seems apparent 
that the use of a Fourier analysis of a meteorological field presents a 
very valuable tool for atmospheric research. With the advent of the 
electronic computer and machine -analyzed contour charts, the use of har- 
monic analysis is an effective and economically feasible approach to the 
computation of synoptic aids for the practising meteorologist. 
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